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The continuum limit and scaling properties of an asymptotically free field theory regularized on a random lattice 
are compared with those on a regular square lattice. We work on random lattices parametrized by a degree of 
"randomness" k. We show that the continuum limit exists and diff'erent k are related by a finite renormalization. 



1. INTRODUCTION 

The lattice regularization was introduced to 
study the non-perturbative aspects of field the- 
ories [1]. To extract continuum physics informa- 
tion from the numerical Monte Carlo, the simu- 
lations must be performed in a region where the 
correlation length ^, the lattice spacing a and the 
lattice size L satisfy 1 ^ ^/a <^ L. How strong 
these inequalities are depends on the theory we 
are studying and the particular lattice used. The 
random lattice was put forward to weaken the 
above inequalities [2]. Indeed, the isotropic pro- 
perties of the random lattice should expedite the 
continuum limit. 

However, the random lattice has been little 
exploited for numerical studies of field theories. 
On the other hand, no much work has been done 
to check whether its continuum limit is the same 
than that of the usual regular square lattice. In 
this work we have considered the 0(3) non-linear 
(T-model in two dimensions in order to address the 
above questions. 

To construct the random lattice, we followed 
the procedure of T. D. Lee et al. [2,3]. The only 
new ingredient is the introduction of a degree of 
"randomness" k [4]. The sites of the lattice are 
the centers of hard spheres, the radius of which is 
a/{2K). These hard spheres are randomly located 
and their relative distance follows the Poisson di- 



stribution. At small values of k the lattice turns 
out to be locally less random, [i.e.: link lengths, 
plaquette areas, coordination numbers, etc. have 
smaller variance) [4] . 

2. THE 0(3) cr-MODEL ON THE RAN- 
DOM LATTICE 

Following [5] , the action of the model on a ran- 
dom lattice can be written as 



'^"^'"' = ^E^«H'^"'«-'^"'^)'- 



(1) 
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In this expression, g is the coupling constant. The 
matrix \ij is the weight that allows Eq.(l) to have 
the correct naive continuum limit. It is defined as 
the ratio between the lengths of the dual link and 
the link joining the sites i and j. 4>i is the value 
of the field 4> at the site i. In a recent paper, [4], 
it has been shown that this model on a random 
lattice is asymptotically free and the renormali- 
zation group invariant scale A depends on k, the 
"randomness" degree parameter. 

Beside these perturbative results, an explicit 
numerical simulation gives us an insight about 
non-perturbative features such as the mass gap 
and the topological content, which can be com- 
pared to the continuum results [6-8]. 

To measure the mass gap, we computed the 
2-point correlation function at zero spatial mo- 
mentum and fitted the result with an hyperbolic 



As a definition of the topological charge opera- 
tor on the random lattice we propose 

Qi = -^ fa6c / ^ 'i't <i'j <l>k ■^ij ■^ik '^ijk, (2) 

A ATT in: < ^ ■* 



continuum limit is reached by fitting the exponent 
r\ in the 2-loop renormalization group prediction 



327r uji 



iyk 



where Qi is the topological charge density at the 
site i, Wi the volume of the Voronoi cell around 
this site [2] and Gij}. is defined as 



Tijh = ihj + dijj X Uik + dikj 



(3) 



In this expression, lij is the link vector joining 
the sites i and j and dij is the vector joining the 
center of the above link with the center of the as- 
sociated dual link [5] . Finally, z is the unit vector 
orthonormal to the plane of the lattice, oriented 
as z = X X y. One can prove that Qi has the cor- 
rect naive continuum limit. The total topological 
charge Q of a configuration is obtained simply by 
Q = ^^ Qi. On a regular lattice, our topological 
charge coincides with the definition of reference 
[9]. 

3. NUMERICAL RESULTS 

In the Monte Carlo simulation we used the 
Wolff updating algorithm [10]. Depending on 
/3 = 1/g, 2-^lOx 10'* initial sweeps are suffi- 
cient to thermalize the configuration. In order to 
compare the performance on each lattice, we repe- 
ated the same measurement on a regular square 
lattice and on a k = 100 and k = 1.3 random 
lattices. These values of k correspond to rather 
opposite levels of "randomness", the first being 
completely irregular and the second almost regu- 
lar. We did not average the result among several 
random lattices since we have used a large enough 
number of sites. On the regular lattice we used 
the standard action, and the action of Eq.(l) on 
the random lattice. Here, we report the results 
obtained with lattices ranging from 200^ to 400^ 
sites. 

3.1. Mass measurement 

For every value of /3, 12000 updatings were 
made. On each of them the Wolff's improved 
estimator was used to compute the 2-point corre- 
lation function [11]. We checked that the correct 



Monte Carlo data = — (0 e"'"^. 
A 



(4) 



The universal beta function dictates r] = ( = 27r. 
We found that r] is equal to 27r within 1% on re- 
gular and random lattices independently of the k 
used. However, we could not measure ( which is 
very sensitive to the non-universal terms. Cor- 
recting the data for finite size effects as shown in 
[12,13] did not change this situation. 

In figure 1 we report the ratio between A pa- 
rameters for both K = 100 and k = 1.3. The 
average values are 1.29(8) and 0.86(6) respecti- 
vely, and do not change when increasing the lat- 
tice size. The theoretical calculated values are 
1.8(2) and 0.92(2) respectively [4]. We think that 
the disagreement between the calculated and me- 
asured ratios for k = 100 may be due to large non- 
universal terms in the beta function. On the other 
hand, we see a rather good agreement in this ra- 
tio for K = 1.3, suggesting a k dependence of the 
non-universal terms in the beta function. Preli- 
minary results for the 0(4) and 0(8) cr-models 
show a much better agreement with the pertur- 
bation theory predictions [14]. 




Figure 1. Ratio of A parameters for k = 100 
(dashed line and circles) and k = 1.3 (solid line 
and squares) calculated on a 200^ lattice. 



3.2. Topological Susceptibility 

We also measured the topological susceptibi- 
lity, defined as 



X 



V 



{Q'), 



(5) 



where V is the lattice volume, using the cooling 
method adapted to the random lattice action. 
The topological charge was measured after 30 co- 
oling steps for each uncorrelated configuration [9] . 
It is interesting to notice that the same clustering 
towards integer values observed after the cooling 
on regular lattices holds here too. We checked 
that this measurement was unchanged within er- 
rors after 40 and 50 cooling steps. We also verified 
that the result for the topological susceptibility is 
insensitive to the rounding of the previous value 
of the topological charge to the nearest integer. 

We measured the topological charge on 1000 
uncorrelated configurations separated by 500 up- 
datings and checked that there was no autocor- 
relation at all. The Monte Carlo data provided a 
value for the r] parameter in the exponent of the 
renormalization group prediction a ~ 20% away 
from the perturbative value 27r for both random 
and regular lattices. This discrepancy diminishes 
when using larger lattices. 

An analysis, analogous to that of figure 1, ap- 
plied on the topological susceptibility data, pro- 
duces similar results: 1.37(3) and 0.94(2) for 
K = 100 and k = 1.3 respectively. 

4. CONCLUSIONS 

By varying the k parameter, we can construct 
random lattices with different levels of "random- 
ness". We checked that they correspond to dif- 
ferent regularizations so that any average among 
random lattices must take into account this fact. 
Our Monte Carlo data support the scenario of a 
common continuum limit as was put forward in 
[4]. The theory regularized on a random lattice 
presents a mass gap which scales as it should. Mo- 
reover, the Monte Carlo data indicate that there 
is a topological content as well. 

We remark that the non-universal terms in the 
beta function are large enough to prevent the on- 
set of asymptotic scaling (we could not determine 



the coefficient (). In that sense, the random lat- 
tice, compared with the regular lattice, does not 
improve the asymptotic scaling. 

A similar analysis for QCD is in progress. 

There, the Arandom is much less than h-regular 

[15]. This fact shifts the scaling window to lar- 
ger betas in a much more effective way than for 
the 0(3) (T-model, possibly allowing a better con- 
trol on the perturbative expansions that mask the 
non-perturbative signal. 
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